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SATELLITE -BASED POSITIONING RECEIVER WITH CORRECTION OF 
CROSS-CORRELATION ERRORS 

Satellite positioning systems employ, for pinpointing, 
5 several satellites that transmit their positions via 
radio signals and a receiver placed at the position to 
be pinpointed, estimating the distances, called pseudo- 
distances, that separate it from the satellites on the 
basis of the propagation times of the satellite signals 

10 picked up and performing the pinpointing operation by 
triangulation. The more precisely the satellite 
positions are known by the receiver and the more 
precise the measurements of the pseudo-distances made 
by the receiver, the more precise the pinpointing 

15 obtained. 

The positions of the satellites are determined on the 
basis of a network of ground tracking stations 
independent of the positioning receivers. These 

20 positions are communicated to the positioning receivers 
via the satellites themselves, by data transmission. 
The pseudo-distances are deduced by the positioning 
receivers from the apparent delays exhibited by the 
received signals relative to the clocks of the 

25 satellites, which are all synchronous. 

Although the precision in knowing the positions of the 
satellites of the positioning system is independent of 
the performance of a positioning receiver, this is not 
30 the case for the precision of the pseudo-distance 
measurements, which depends on the precision of the 
measurements of the signal propagation times at the 
receiver . 

35 Radio signals transmitted by satellites travel large 
distances and, since they are transmitted at limited 
power levels, reach the receivers with very low power 
levels that are buried in radio noise due to the 
activity of terrestrial transmitters which are often 
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powerful and much closer to the receivers than the 
satellites that they must receive. To make it easier to 
receive them, it has been attempted to make them the 
least sensitive possible to narrow-band interference, 
5 by increasing their bandwidths by means of the band 
spreading technique. The current systems, and those 
planned for the near future, for satellite positioning 
use, for the radio signals transmitted by their 
satellites, the technique of band spreading by 

10 modulation with the aid of pseudorandom binary 
sequences, a technique known as DSSS (Direct Sequence 
Spread Spectrum) . This DSSS modulation consists, after 
having arranged the information to be transmitted in 
the form of a sequence of binary elements with a 

15 regular data rate, in multiplying each binary 
information element by a pseudorandom binary sequence 
of markedly faster data rate. The band spreading 
obtained is proportional to the ratio of the data rate 
of the sequence of binary data elements to the data 

20 rate of the pseudorandom binary spreading sequence. 

The information to be transmitted from the satellites, 
once placed in the form of a frequency-spread sequence 
of binary data items by DSSS modulation, are transposed 

25 in the transmission frequency range by modulation with 
a transmission carrier. To make it easier to measure 
the signal propagation times at a positioning receiver 
and to avoid the presence of isolated lines in the 
spectra of the signals transmitted by the satellites, 

30 each pseudorandom binary sequence used for frequency 
spreading consists of binary elements of the same 
duration, taken to be equal to integer multiples of the 
periods of the transmission carriers, whereas the 
various data rates and frequencies used within the 

35 satellites are synchronized and derive from a very 
precise common clock. 

Upon reception, the binary information contained in a 
radio signal from a satellite of a positioning system 



is extracted by two demodulations carried out in an 
intertwined manner, a first demodulation using a 
carrier generated locally by an oscillator controlled 
by a PLL (Phase Lock Loop) for transposing the signal 
5 received into baseband and a second demodulation using 
pseudorandom binary sequences generated locally by a 
pseudorandom binary sequence generator controlled by a 
so-called DLL (Delay Lock Loop) for despreading the 
binary information string present in the received 
10 signal. 

The propagation times of the received signals are 
manifested, at reception, by delays that affect the 
pseudorandom binary sequences present in the received 
15 signals and the carrier modulating the received signal. 

The delays affecting the pseudorandom binary sequences 
are accessible, modulo the duration of one of their 
binary elements, at the level of the feedback control 

2 0 signals of the DLLs. The delays observed by these loops 
allow unambiguous measurements, or those of low 
ambiguity, of the propagation times of the pseudorandom 
binary sequences since the numbers of complete 
pseudorandom sequences passing during the signal 

25 journeys is relatively small. This is referred to as 
code measurements . 

For example, in the case of the GPS (Global Positioning 
System) satellite positioning system, the shortest 

30 pseudorandom binary sequence, that used for satellite 
signal spreading of the C/A (Coarse/Acquisition Code or 
Clear/Acquisition Code) type, is composed of 1023 
binary elements with a data rate of 1023 MHz and a 
duration of one millisecond. Its total duration 

35 corresponds to a journey of 300 km for a radio wave and 
allows measurements of distance modulo 300 km. The 1 
microsecond duration of each of its binary elements 
permits a precision of the order of 0.1 microseconds in 
the measurement of its delay at reception, 
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corresponding to a 30 meter journey in the case of a 
radio wave. The ambiguity in the pseudo-distance 
measurements obtained from the pseudorandom binary 
sequence of a C/A code, due to the fact that 
5 measurements modulo 300 km are involved, is easy to 
resolve as soon as the receiver receives more than four 
satellites, as it is then possible to take various 
points on the same position from different sets of four 
satellites and to retain only the common solution. In 

10 the absence of such a possibility, the ambiguity may 
also be resolved using very rough prior knowledge of 
the position. Such a measurement ambiguity does not 
arise with P-type satellite signals of the GPS system, 
which use, for spreading them, a pseudorandom binary 

15 sequence of 266.41 days' duration, but these signals 
are not freely available to users . 

Figure 1 shows the schematic of a satellite-based 
positioning receiver of the prior art. The receiver 
2 0 comprises a correlator channel 10 driven by the signal 
received, originating on the one hand from the 
positioning satellites visible to the antenna of the 
receiver, and on the other hand from a disturbing 
source . 

25 

The correlator channel 10 comprises a correlation 
path 12 for in-phase and quadrature correlation between 
the signal received Sr and two respective local 
carriers Fi,Fq. These quadrature local carriers (sine, 
30 cosine) are generated by an oscillator with digital 
control of carrier 14 (NCO p) of the receiver. 

The signals I, Q output by the carrier correlation path 
are thereafter correlated in a code correlation path 16 
35 with the local codes, punctual and delta, of the 
satellite considered, provided by a digital generator 
of local codes 18. 

The code correlations are thereafter integrated by a 
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respective integrator 20 so as to provide signals 
IpjIzwQpaQa at the output of the correlator channel 10. 

A carrier loop 22 and a code loop 24, of the receiver 
5 of the state of the art, embodied in a known manner, 
respectively provide on the basis of the signals 
Ip/l^Qp/Qa at the output of the correlator channel 10, 
a carrier speed signal Vpc for controlling the 
oscillator with digital control of carrier (NCO p) 14 
10 generating the two local carriers F lf FQ and a code speed 
signal Vcc for controlling the oscillator with digital 
control of code 18 generating the local code, punctual 
and delta, for the code correlation path. 

15 Satellite-based radio navigation requires means on the 
ground (grourtd segment) in order to monitor and correct 
the signals emitted by the satellites. The ground 
stations in particular use means of reception that 
provide code and carrier measurements. These 

20 measurements must be extremely accurate since they 
contribute to the ultimate performance of the system. 

These measurements are marred by errors due to cross- 
correlations between the signals emitted by all the 

25 satellites visible from the ground stations. 
Specifically, although the spreading codes used are 
designed to distinguish the signals of the satellites 
by correlation, their decorrelations are not perfect on 
account of their limited length and of the Doppler due 

30 to the motion of the satellites. These errors may be 
perturbing when the codes are of short period (C/A 
codes) and when the speeds between the ground and the 
satellites are low. This is the case in particular for 
the GEO and SBAS (Space Based Augmented Systems, type 

35 WAAS or EGNOS) satellites: the small variation in the 
Doppler implies that these errors become slowly varying 
biases that cannot be filtered. They may reach several 
meters . 



The idea of the invention is based on the use of 
additional correlation channels over and above the 
correlation channel of the signal received from a 
satellite so as to estimate in real time the cross- 
5 - correlation errors, code-wise and carrier-wise, between 
the satellite concerned and any other satellite; which 
we track moreover on other channels and the position of 
the code and phase of whose carrier we therefore also 
know. These estimated errors may thus be corrected very 
10 simply in the tracking loops. 

To this end the invention proposes a satellite-based 
positioning receiver with correction of inter-satellite 
cross-correlation errors, the receiver comprising a 
15 correlation channel Cii of order i per satellite 
received, with i=l,2,..N, N being the number of 
satellites received, each correlator channel Cii 
having: 

a carrier correlation path, in-phase and 
20 quadrature, between the signal received Sr and two 
respective local quadrature carriers (sine, cosine) 
generated by an oscillator with digital control of 
carrier (NCO p) ; 

a code correlation path based on the signals I, 
25 Q output by the in-phase and quadrature carrier 
correlation path, with the local codes of the satellite 
received, provided by a digital generator of local 
codes ; 

an integrator for providing, for each local 
30 code, signals I c Q c at the output of the correlator 
channel Cii of the satellite received, c designating 
each of the local codes, 

characterized in that it comprises, for each correlator 
channel Cii of the satellite received as many 
35 additional correlator channels Cix as additional 
satellites received with x=l,2,...N and x different 
from i, and in that the local codes of the satellite 
received are correlated with the local codes of the 
other additional satellites Cix. 



In one embodiment of the receiver according to the 
invention the local codes of the satellite received for 
the code correlation path are a punctual and delta 
5 code. The code correlation path in fact comprises two 
correlation paths: 

a punctual path (I p ,Qp), 
a delta path (I a ,Qa). 

In another embodiment, the local codes of the satellite 
received for the code correlation path are an early, 
punctual and delta code. The code correlation path in 
fact comprises three correlation paths: 
- an early path (I A ,Q A ), 

a punctual path (Ip/Qp), and 

a late path (Ir,Qr), the delta path being 
reconstituted from the early path minus the late path 
by the formulae: 
1a = 1a " Ir 
Qa = Qa " Qr 

In a first embodiment of the receiver according to the 
invention, receiving the N satellites, the receiver 
comprises N reception subsets Si. Each subset Si of 
2 5 rank i comprises the correlator channel Cii of the 
signal of the satellite received of order i, with i = 
1,2,...N, and N-l additional correlator channels 
Cil, . .Cix, . . .CiN for the additional satellites received 
with x = 1,2,... N and x different from i. 

30 

Each received signal correlator channel Cii is driven 
by its reception input Er by the signal received Sr. 
Each of the additional correlator channels of a subset 
Si receives respectively, on the one hand, at its 
35 received-signal input Er, a local signal Slox resulting 
from the modulation of the local carrier F lx by the 
punctual local code Cpx of the correlator channel Cxx 
of the satellite received of order x, and on the other 
hand, at its local carrier and local codes local 
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inputs, the respective local quadrature carriers F^Fqi 
and the local codes Cpi and Ai of the correlator 
channel Cii of the signal received from the satellite 
of order i. 

5 

The invention will be better understood with the aid of 
an exemplary embodiment of a receiver according to the 
invention receiving more than two satellites with 
reference to the appended drawings, in which: 

10 

figure 1, already described, represents a 
receiver of the state of the art having a correlator 
channel; 

figure 2 shows a subset of a receiver according 
15 to the invention receiving N satellites; 

figure 3 shows a receiver, according to the 
invention, for three satellites; 

figure 4 shows a correlators channel operating 
in baseband; 

20 - figure 5 shows a subset of the receiver 

according to the invention operating in baseband; 

figure 6 shows a baseband receiver, according 
to the invention, for three satellites. 

25 Figure 2 shows a subset of a receiver according to the 
invention receiving N satellites. 

The receiver comprises N reception subsets for the N 
satellites received. Each subset Si of rank i, with 

30 i=2,3,...N, comprises a correlator channel Cii for a 
satellite received Sati of order i and N-l additional 
correlators Cii, Cix, . . . CiN for the additional 
satellites Satl, . . . Satx, . . . SatN, with x different from 
i. Each of these additional correlator channels Cix 

35 receives respectively, on the one hand, at its 
received-signal input, a local signal Slox resulting 
from the modulation of the local carrier F ix by the 
punctual local code Cpx of the correlator channel Cxx 
of the signal received from the satellite of order x, 



and on the other hand, at its local-carrier and 
local-codes inputs, the respective local quadrature 
carriers Fn, F Q i and the local codes, punctual Cpi and 
delta Ai, of the correlator channel Cii of the signal 
5 received from the satellite of order i. 

Like the received-signal correlator channel Cii, each 

additional correlator channel of rank x in the subset 

Si, if x=l,2,...N comprises: 
10 - the in-phase and quadrature carrier correlation 

path 12 between the signal received and two respective 

quadrature local carriers (sine, cosine) ; 

the code correlation path 16 based on the 

signals I, Q at the output of the in-phase and 
15 quadrature , carrier correlation path with the punctual 

Cpi and delta Ai local codes of the satellite of order 

i; 

an integrator for providing signals Ipi X , 
Iaix/Qpix/ Qmx at the output of the correlator channel. 

20 The integrator of the received signal correlator 
channel Cii provides the signals Ipn, Iml, Qpii, Qaii. 
The subset Si furthermore comprises: 

an oscillator with digital control of carrier 
OPi (NCO p) for providing local carriers Fn, F Qi for the 

25 N correlators of the subset Si considered and a digital 
generator of local codes OCi (NCO c) for providing the 
local codes, punctual Cpi and delta Ai, for the N 
correlators of the subset Si considered; 

a multiplier Mi providing, for the other 

30 subsets Sx of the receiver, a local signal Sloi, 
resulting from the modulation of the local carrier Fn 
by the punctual code Cpi of the subset considered Si, 
so as to perform the correlation of the code modulated 
by the carrier of the satellite considered with the 

35 codes modulated by the carrier of the other satellites; 

a correlation corrector CRi providing on the 
basis of the signals I Pix , I Aix , Q P i X , Qaix (x taking, for 
these signals I P i X , I A ix/ Qpi X/ Qai*/ the values 1 to N) at 
the output of the N correlator channels of the subset 
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considered Si, and signals I Pxx , Q pxx output by the 
received-signal correlator channels of the other 
subsets Sx, with x different from i, corrected signals 
Ipi' ,Im' 9 Qpi' QAm' - 
5 a carrier discriminator DPi providing through a 

carrier loop corrector CBPi a control signal Vcpi for 
the oscillator with digital control of- carrier (NCO p) 
so as to provide local carriers F Ii# F Q i for the N 
correlators of the subset Si considered; 
10 - ^a code loop discriminator DCi providing through 

a code loop corrector CBCi a control signal Vcci for 
the digital generator of local code OCi (NCO c) for 
providing the local codes, punctual Cpi and delta Ai 
for the N correlators of the subset Si considered. 

15 

We shall deal hereinbelow, by way of example, with the 
case of a receiver according to the invention, 
configured to receive three satellites (N=3) . 

20 Figure 3 shows a receiver of three satellites 
comprising a first SI, a second S2 and a third S3 
reception subsets having three correlator channels 
each. The reception subsets SI, S2 and S3 comprise the 
same elements as the detailed subset of figure 2. 

25 

The first SI, second S2, and third S3 subsets of the 
receiver of figure 3 respectively comprise a first Cll, 
a second C22 and a third C33 signal correlator channels 
driven at their reception input Er by the signal Sr 
30 received by the receiver, each subset furthermore 
comprising : 

in the first subset SI, two other additional 
correlator channels C12 and C13 driven respectively at 
their reception input by local signals Slo2, Slo3 
35 emanating respectively from a multiplier M2 and from a 
multiplier M3, the signal Slo2 resulting from the 
modulation of the local carrier Fi 2 by the punctual code 
Cp 2 of the second satellite and the signal Slo3 
resulting from the modulation of the local carrier F r3 
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by the punctual code Cp 3 of the third satellite; 

in the second subset S2, two other additional 
correlator channels C21 and C23 driven respectively at 
their reception input by local signals Slol, Slo3 
5 emanating respectively from a multiplier Ml and from 
the multiplier M3, the signal Slol resulting from the 
modulation of the local carrier Fn by the punctual code 
Cpi of the first satellite; 

in the third subset S3, two other additional 
10 correlator channels C31 and C32 driven at their 
reception input by the local signals Slol, Slo2 
emanating respectively from the multipliers Ml and M2 . 

Each correlator of each of the subsets Si comprises: 

15 - the in-phase and quadrature carrier correlation 

path 12 (see detail of the correlator in figure 2) 
between the signal at their reception input and two 
respective quadrature local carriers (sine, cosine), 
Fii,F Q i for the first subset SI, Fi 2 ,F Q2 for the second S2 

20 and Fi 3 ,F Q3 for the third S3, these carriers being 
generated respectively, for each of the subsets S1,S2, 
and S3 by a first OP1, a second OP2 and a third OP3 
oscillators with digital control of carrier (NCO p) ; 

the code correlation path 16 based on the 

25 signals I, Q at the output of the in-phase and 
quadrature carrier correlation path with the local 
codes, punctual Cpl,Cp2,Cp3 and delta Al,A2,A3 of the 
satellites respectively of order 1, 2, 3, provided by a 
digital generator of local codes OC1, OC2 and OC3 

30 respectively for each subset; 

an integrator per correlator channel for 
respectively providing signals I P i x , I A i x , Q P i x , Qaix at the 
output of the correlator channel Clx; Ip2x> Ia2x* Qp2x/ Qa2x 
at the output of the correlator channel C2x and 

35 I P3x , I A3x , Q P3x , Q A3x at the output of the correlator channel 
C3x, with x=l,2,3. 

Each subset of three correlators comprises: 

a corrector Crl,Cr2,Cr3 of correlations 
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providing on the basis of the signals Ipi X , I^ix^ Qpix> Qaixj 
with i=l,2,3 (and x=l,2,3), at the output of the 3 
correlator channels of the subset considered S1,S2,S3 
and of the signals Ipxx/Qpxx/ at the output of the 
5 received-signal correlator channels (of order x) of the 
other subsets Sx, of the corrected signals 
Ipi' # Iai' f Qpi' ,Qai' at the output of the first corrector 
Crl, I P2 ' , Ia2' f Qp2' , Qa2' at the output of the second 
corrector Cr2, I P3 ' , I A3 ' , Q P3 ' , Q A3 ' at the output of the 

10 third corrector Cr3, the signals Ipxx/Qpxx at the output 
of the received-signal correlator channels, driving the 
correctors, being the signals I P 22, Ip33# Qp22# Qp33 for the 
corrector Crl, Ipn, Ip33/ Qpiir Qp33 for the corrector Cr2 and 
Ipn, Ip22aQpii/QP22 for the corrector Cr3, 

15 - a carrier discriminator DP1,DP2,DP3 

respectively providing through a carrier loop corrector 
CBP1,CBP2,CBP3 a control signal Vcpl, Vcp2, Vcp3 for the 
respective oscillator with digital control of carrier 
OPl,OP2,OP3 (NCO p) so as to provide local carriers 

20 Fn,F Q i, for the first subset SI, Fi 2 ,Fq 2 , for the second 
subset S2 and F I3 ,F Q3 for the third subset S3; 

- a code loop discriminator DC1,DC2,DC3 
respectively providing through a code loop corrector 
CBC1,CBC2,CBC3 a respective control signal 

25 Vccl, Vcc2 , Vcc3 for the digital generator of local codes 
OCl,OC2,OC3 (NCO c) so as to provide the local codes, 
punctual and delta, Cpl,Al for the three correlators of 
the first subset SI, Cp2,A2 for the three correlators 
of the second subset S2 and CP3,A3 for the three 

30 correlators of the third subset S3. 

The receiver of figure 3 is configured to perform the 
following corrections : 

35 For the satellite 1: 

On the punctual path: 
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IpV - Ip11 " IP22 • IP12 • 2/T - Ip33 ■ ip13 ♦ 2/T 
Qpi' - Qpn - Ip22 . QP12 • 2/T - Ip33 ■ Qpi3 • 2/T 



on the delta path: 

Ul' = Ul1 - IP22 • Ul2 2/T - IP33 , Ia13 . 2/T 
5 Q A 1* = Q A 11 - IP22 - Qa12 2/T - Ip33 . Q A 1 3 . 2/T 

i.e. in complex notation, with j 2 =-l: 

IP1*+ jQpi - IP11+ JQph - IP22 (lp«+ JQP12). 2/T- lp 33 (Ipi3 +jQpia). 2/T 
Lm' + JQai'= JQah - Ip22 (Ui2 + ]Qai 2 ). 2/T- Ip33 (IA13 +jQAl3). 2/T 

with —= foocal signal (t)) 2 dt,T duration of integration of 
2 o 

10 the integrator 20 in figure 1. 



Proof : 



Remark: in the figures and, with the aim of making it 
15 simpler to read them, we shall write: 



Punctual local code sa teiiite i = Slol 
Punctual local code sa teiiite 2 = Slo2 
Punctual local code sa teiiite 3 = Slo3 

20 

By construction: 
T 

Qp { = J[ *gnal received® ]\localcode ^//^C). Quadrature local carrier ^^^(01* 

i.e., by expanding and neglecting noise: 

25 

7 

I Pl = |[/4,xo^ 1 0 + r i ).sin(^/ + ^ 1 ) + /4 2 ^e 2 (/ + rJ.sin(^+^ 2 ) + 

0 

A 2 .code 2 {t + r 2 ). sin( cot + <p % )][ code , (t + r, ). sin( <2>f + <p, )] r# 

7" 

Qp\ = jMi ( / + r i )• sin(firf + ) + A l .code 2 (/ + r 2 ). sin(tftf + ^ 2 ) + 

0 

A 2 .code 2 (t + r 2 ). sin(o)t + p 3 )] [code, (r + r, ). cos(« + ^ )] rf/ 
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t 

./,,, = J(^ r cWe l (/ + r l ).sin(aj/ + ^ ) )][co^,(/+rJ.sin(a>/ + ^ l )]^ 

o 

T 

- A 2 j[code 2 (t + r 2 ) .sin(tf>r + #> 2 )][corfe 1 (f + r, ). sin(firt + )] dt 

o 

- A> j[co<&? 3 (/ 4- r 3 ).sin(<atf + (p 2 )][code 1 (t + r, ).sin(#/ + )] 



2pi = J[^i code y (t + ). sin(atf + )][ce> Je, (t + Tj ).cos(atf + <p x )] </f 
o 

7" 

- A 2 j[code 2 (t + r 2 ).s'm(cot + <p 2 ))[code x (t + Tj ).cos(6tf + <p } )] A 

0 

r 

- j[code 2 (t + r 3 ). sin(<y/ + p 3 )][code } (t + r, ).cos(tf>/ + <p x )] eft 



0 

We would ideally like: 



T 
r 

°« ' = J I*"™ BuMtoflMW"* satellite!^ " retoM/ ^W//ifd W] * 



5 i.e. by expanding; 



f /'i ' - J[ ^ i ' corfe i (' + r > )• sin(<y/ + )] [code } {t + r, ). sin(<ar + )] rf/ 

o 

T 

Qr\ ' = J[- 4 i (< + T » ).sin(aar + <p, )][corfe, (r + r, ). cos(<uf + )] dt 



/,.,* = /,,, - A 2 j[code 7 (/ + r 2 ). s\n(o)t + <p z )][code x {t + r, ). sin(o;/ + <p { )] dt - 

o 

r 



^3 J[«wfe 3 (f + r 3 ). sin( cot + p 3 )][ cdwfe , (/ + r t ). sin( cot + )] dt 
o 

r 

£?/>i ' = Qpi~ A i \[code t (t + t % ).sm(cot + <p 2 )][code } (/ + r, ). cos(d* + <p x )] dt - 

o 

T 

/ij J[co<fe 3 (/ + r 3 ).sin(6X + p 3 )][c0*fe, (/ + r x )xos(cot + q> x )] rff 

0 

given that: 
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r 

I m = j[ctfde 2 (/ + r 2 ).sm(^/ + ^ 
o 
r 

(?/>i2 = \i c °de 2 (t + r 3 ). sin(a>r + p 2 )][co^ (r + r, ).cos(*tf + ^, )] rfr 

0 
7' 

0 

r 

£?/mj 55 Jfcorfe, (/ + r a ). sin(M + ^ )][corfe 1 (f + r, ), cos(o» + )] A 

0 

and, neglecting the inter-satellite cross-correlation 
terms, and assuming that the local carriers are in 
phase with the carriers received: 

I m = j[A z .code 2 (t + r 2 ). sin(#/ + )] [cotf<? 2 (t + r 2 ). sin(atf + ^ 2 )] (ft « — ^4 2 



= j[ ^ 3 .corfe, (* + r 3 ). sin(<y/ + p 3 )][c0<fe 3 (/ + r 3 ). sin(ot + #> 3 )] = — 

We do indeed obtain the formulae proposed when we 
replace the correction terms with the corresponding 
terms I P i X and Q P i X . 

10 

Likewise : 

r 

= -A 2 J[cot/e 2 (/ + r 2 ),sin(^ + ^ 2 )][^/^,(/ + rJ.sin(6>/ + ^ I )]^- 

0 

T 

A y j[code 3 (t + r 3 ). sin(a# + <p 3 )][delta x (t + r, ). sin(<# + )] dt 
o 

r 

Qai ' - 2ai - ^ 2 JO><fe 2 (/ + r a ). sin(<u/ + p 2 (r + r, ). cos(art + cp x )] <fr - 

0 

r 

4 J[c<?tfe 3 (/ + r 3 ). sin(<a>/ + #> 3 )][delta { (t + r, ). cos(cot + q> x )] A 



r 

0 
T 

Qsm = J[codc 2 (/ + r 2 ).sin(otf + ^ 2 )][delta x {t + r, ), cos(<a/ + )] rf/ 

0 

15 
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In the case where the local carriers are not entirely 
in phase with the carriers received, on account of the 
dynamics (carrier, receiver clock, satellite) we show 
that : 

5 for the first satellite Satl: 

- for the punctual path: 

Ipi' - Ipi1 — ( Ip22 > 'P12 " QP22 • QP12)- 2/T — ( Ip33 . Ipi3 - Qp33 , QP13). 2/T 
Qpi* = Qpn - ( Ip22 -QP12 + Qp22 • IP12)- 2/T - ( Ip33 • Qpi3 + Qp33 • Ipi3)- 2/T 

10 - on the delta path: 

Ul' = 'a11 - ( IP22 • Ul2 " QP22 • Qa12). 2/T - ( l P33 . I A1 3 - Qp 3 3 . Q A 13). 2/T 
Qa1* = Qa11 ~ ( 'P22 Qa12 + Qp22 . Iai2)- 2/T - ( | P33 . Q A13 + Qp 33 . I A13 ). 2/T 

i.e. in complex notation, with j 2 =-l: 

15 

I Pi' + j Op\ = Ipii + j Qpil — ^ on x different from i (lpxx + jQpxx)(lpix + jQpix)2/T 
Iai' + j Q^' " I/Mi + ] QaII - 2 on x different from! (Ipxx + jQpxx)OAlx+ jQAix)2/T 

for the second satellite Sat2 : 

| P2 ' + jQp2 = Ip22+ jQp22-0P11 + jQpn)(IP21+ jQp2l)/T-(lp33 + J'QP33)(IP23 + jQP23)/T 
l A2 ' + jQ*2' = IA22+ jQ A 22-(lpi1+ jQpil)(U21+ jQA2l)/T-{lp33+ jC>P33)(lA23+ ]Q^)IT 

20 

for the third satellite Sat3: 

IP3 + jQp3 - IP33+ jQp33"(lpi1 + jQpil)(lp31 + jQp3l)/T-(I P2 2- jQp22)(lp32+ jQp32)/T 
Ia3' + jQA3 - U33+ jQA33 "(Ipi1 + JQpu)(Ia31 + jQA3l)/T -(Ip22~ jQ P2 2)(lA32+ jQ^IT 

25 

Generalization: 
On the punctual path: 
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I Pi = Ipii 53 on x different from i ( Ipxx • tax ~ Qpxx • Qpix)- 2/T 
Qpj' ^ Qpit — E on x different from i ( Ipxx -Qpix + Qpxx • Ipix)- 2/T 

On the delta path: 

L> = I Ail ~ S on x different from i ( Ipxx « Uix " Qpxx * Qaix)* 2/T 
Qa! 23 QaSi — £ on x different from i ( Ipxx *Qaix + Qpxx • Uix)* 2/T 

5 

i.e. in complex notation, with j 2 =-l: 

Ipi' + j Qp\ 3 Ipil + j Qpll ~ 2 on x different from i (lp xx + jQpxx)(lpix + jQpix)2/T 
Iai' + j CW = lAIi + j QaII - 2 onxdfferentfromi (| Pxx + jQpxx)OAlx+ jQAix)2/T 

10 

In order to make the notation for the indices more 
systematic, the index ii labels the correlator channel 
Cii of the subset Si which handles the signal received, 
different from the other correlator channels Cix of the 
15 'subset Si which, for their part, handle the local 
signals from the other satellites of respective order 
x, emanating from the correlator channels Cxx of the 
other subsets Sx. 

20 In a variant of the receiver with correction of cross- 
correlation errors, according to the invention, the 
correlator channels are driven in baseband with signals 
I and Q. 

25 Figure 4 shows a correlator channel 50 operating with a 
signal received Br in baseband. As in the case of the 
correlator channel 10 of figure 1, the baseband 
correlator channel 50 comprises an in-phase and 
quadrature correlation path 52 between the baseband 

30 signal received, in the form of two signals I and Q in 
quadrature, and two respective local carriers Fi , Fq . 
These local quadrature carriers (sine, cosine) are 
generated by an oscillator with digital control of 
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carrier 54 (NCO p) of the receiver. 

The signals I,Q output by the carrier correlation path 
are thereafter correlated in a code correlation path 56 
5 with the local codes, punctual Cp and delta A, provided 
by a digital generator of local codes 58. 

The code correlations are thereafter integrated by a 
respective integrator 60 so as to provide signals 
10 Ip/Ia,Qp/Qa at the output of the correlator channel 50. 

Figure 5 shows a subset of rank i of the receiver 
according to the invention operating in baseband. 

15 As in the receiver of figure 2 operating with the 
signals received, the baseband receiver comprises N 
reception subsets for N satellites received. Each 
subset Si of rank i, with i=2 , 3, . . . N, comprises a 
correlator channel Cii for a satellite received Sati 

20 and N-l additional correlators Cii, Cix, . . . CiN for the 
additional satellites Satl, Satx, . . . SatN, with x 
different from i. The correlator channel Cii and the 
additional channels having the structure of the 
baseband correlator channel of Figure 4. The subset Si 

25 furthermore comprises: 

an oscillator with digital control of carrier 
OPi (NCO p) for providing local carriers Fu^Fqi for the 
N correlators of the subset Si considered and a digital 
generator of local codes OCi (NCO c) for providing the 

30 local codes, punctual Cpi and delta Ai, for the N 
correlators of the subset Si considered; 

a first Mli and a second MQi multipliers 
providing for the other subsets of the receiver a first 
SLIi and a second SLQi local signals resulting from the 

35 modulation of the quadrature signals Fn and F Q i of the 
local carrier by the punctual code Cpi of the subset 
considered, so as to perform the correlation of the 
code modulated by the carrier of the satellite 
considered with the codes modulated by the carrier of 
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the other satellites; 

a correlation corrector CRi providing on the 
basis of the signals Ipj. x/ Iaix, Qpix,Qaix at the output of 
the N correlator channels of the subset considered Si 
5 and signals I PX x, Qpxx at the output of the 
received-signal correlator channels of the other 
subsets Sx, with x different from i, corrected signals 

I Pi' t lAi' i Qpi/ QAi' • 

a carrier discriminator Dpi providing through a 
10 carrier loop corrector CBPi a control signal Vcpi for 
the oscillator with digital control of carrier (NCO p) 
so as to provide local carriers Fii,F Q i for the N 
correlators of the subset Si considered; 

a code loop discriminator DCi providing through 
15 a code loop corrector CBCi a control signal Vcci for 
the digital generator of local codes OCi (NCO c) for 
providing the local codes, punctual Cpi and delta Ai 
for the N correlators of the subset Si considered. 

20 Figure 6 shows a baseband receiver for three satellites 
Satl,Sat2 and Sat3 comprising a first SI, a second S2 
and a third S3 reception subsets having three 
correlator channels each. The reception subsets S1,S2 
and S3 comprise the same elements as the detailed 

25 subset of figure 4 operating in baseband. 

The receiver of figure 6 is configured to perform the 
same corrections as those of the receiver of figure 3 
except that T/2 is replaced by T in the correction 
30 formulae. 

For example for satellite 1: 

IP1 ,+ jQpi' = IP11+ jQpil -0p22+ jQp22)(lp12+ jQpi2)/T -(Ip33+ jQp33)(lpi3+ jQpi3)/T 
Lm'+ JQai' = jQ*11 -flP22+ jQp22)(Ul2+ ]Q A 12)/T -(l P3 3+ jQp33)(Ul3+ jQ A 13)/T 

35 

In a variant configuration, the receiver according to 
the invention uses three code correlators: 
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an early path (I A/ Q A ) , 
a punctual path (I P , Q P ) , and 
- a late path (I R , Qr) , the delta path being 
equivalent to the early path minus the late path, 

5 

In another variant of the receiver according to the 
invention, the same method is applied to the delta path 
reconstituted at the output of the correlators by the 
formulae : 

10 laix = lAix~lRix 
Q&ix = QAix~QRix 

The cross-correlations are calculated twice in the 
above description. In fact, it is possible to economize 
15 on the correlators by virtue of the formula below: 

For the first satellite Satl, we calculate 
(Ip, Ia, Op/ Qa) 12 and (I P , Ia,Qp,Qa) 13 in addition to 
( Ipf Ia, Qp, Qa) ii 

20 

lpi'+ JQpi' = Ipi1 + J'QP11 -0p22+ jQp22)(lpi2+ j QP12)/T-(I P 33+ j Qp33)(lpi3+ jQpisVT 
jQAV = 1*11+ ]Q A 11 "(IP22+ jQp22)(l A 12+ jQAl2)/T-(l P3 3 + jQp33)(lA13+ \Oi^)IJ 

For the second satellite Sat2, we calculate 
Up, Ia/Qp/Qa)23 in addition to ( I P , I A , Qp, Qa) 22 

25 

lp 2 '+ jQ P2 ' = Ip22+ jQp22-('P11 + JQP11)(IP12 * jQp12)/T-(lp33+ jQp33)(lp23 + jQP23)/T 
Ia2'+ ]Q& = 1*22+ jQA22 +(IP11+ jQpil)(lA12 " jQAl2)/T-(l P3 3+ jQp33)(U23+ \Q^)fT 



For the third satellite Sat3, nothing is calculated in 
addition to ( I P , Ia, Qp, Qa) 33 

30 

Ip3 + JQP3 = IP33+ jQp33 -(IP11+ jQpil)(lpi3- jQpi3)/T-(lp22+ jQp22)(lp23" jQp23)fi~ 
Ia3' + jQA3' = lA33 + jQA33+(lpi1+ jQpil)(lA13-jQA13)/T+(lp22 + jQp22)(l A 23- jQA23)/T 



Generalizing, for x > i: 
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Ipxi = + Ipix 
Qpxi = ~ Qpix 
^Axi = ~ ^Aix 
Qaxi = + Qaix 

5 

To summarize, instead of having N 2 -N correlator 
channels in addition we have (N 2 -N) /2 of them. 

The terms Ipn and Q Pii in the formulae, estimation of 
10 the complex amplitude of the signals received 
respectively from the satellites i, take no account of 
the corrections. In order to improve the accuracy, they 
could be replaced by I P j/ and Q P j/ in the formulae. In 
this case, they become: 

15 

Ipi' + j Qp\ ~ lp|f + j Qpii — 2 on x different from i (W + jQpx')(lpix + ]Qpix)2/T 
Ur + j CW - Uli + j QaII - 2 on x different from i (lpx' + jQpx)(Ulx + jQAIx)2/T 

The problem which then arises is that the application 
of the formulae demands inputs I P i' and Q P i' which are 

20 themselves the outputs of the calculation. In order to 
overcome this, in a variant of the receiver according 
to the invention, at each iteration of the calculation, 
the corrected terms I P j/ and Q P i' of the previous 
iteration may be used, initializing the calculation 

25 with uncorrected terms I Pii and Q P ii, after the 
acquisition and convergence phase: 

( lp|' + j Qp\ )n = ( Ipii + j Qpil )n - Son x different from i (lpx' + jCW)n-1 -(lpix + jQpix)n - 2/T 
( W + j Qai' )n = { Ui! + j QAii )n ~ Son x different from i (W + jQpx')n-1 .('Aix + jQAix)n - 2/T 

30 The iteration indexed by n may be either over time, 
corresponding in each instance to new data, or that of 
a recursive calculation converging to the ideal 
solution. In a variant of the receiver according to the 
invention, when the signal received is filtered 

35 (limited spectrum), it is possible and recommended to 
apply the same filtering to the local signals. On the 
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other hand, this makes it necessary to install a filter 
for each local signal, hence one per satellite, 
contrary to the signal received, of which there is only 
one . 

5 

We shall see hereinbelow a method of acquisition of the 
signal by the receiver according to the invention. 

A first satellite is acquired, without correction, by a 
10 conventional open-loop search process, well known to 
those skilled in the art. On completion of this process 
we switch to tracking, we deduce therefrom the local 
signal of this first satellite and we correct the 
cross-correlations on the other channels in the search 
15 phase (in open loop) . This makes it possible to acquire 
the weakest satellites (last) while reducing the risk 
of a mistake on account of a correlation with the 
signal from another more powerful satellite. 

20 Each time a new satellite is acquired and tracked, we 
calculate and we apply the cross-correlation 
corrections in respect of the measurements of all the 
other satellites already tracked. 

25 The receiver according to the invention exhibits 
excellent stability. Specifically, because the cross- 
correlation coefficients are appreciably less than 1 
(-24 dB for the C/A codes), the tracking loops are 
stable and converge to a state where there is no longer 

30 any cross-correlation error. 

The receiver according to the invention allows the 
estimation of cross-correlation errors in real time, on 
the punctual and delta aggregate samples I and Q, by 
35 virtue of additional channels, by correlation between 
the local codes of the satellites tracked and the 
correction of the punctual and delta aggregate samples 
I and Q ahead of the code and carrier phase 
discriminators . 
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The receiver according to the invention completely 
eliminates the errors of cross-correlation between all 
the satellites whose signal is tracked, in the steady 
5 state, after a phase of fast convergence. The residual 
errors, due to thermal noise and to the tailing off of 
the loops, depend on the signal-to-noise ratio, on the 
dynamics and on the loop bands. For applications with 
very weak dynamics (ground station) the gain of the 
10 procedure may be very substantial, changing the 
measurement error from a few meters to a few tens of 
centimeters, i.e. a factor of 10. 



the integrator (20) . 



